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NUMERICAL  MODELING  OF  THE 
TEMPERATURE  GRADIENTS  AND  SPECTRA  OF  THE 
PLASMA  ARMATURE  ON  THUNDERBOLT 


INTRODUCTION 


The  puipose  of  this  work  is  to  demonstrate  how,  via  numerical  modeling  of  spectra,  one 
can  infer  some  information  about  the  hot  interior  of  a  dense  plasma  from  the  radiation  which 
is  emitted  principally  by  the  optically  thick  cold  exterior.  The  experimental  data  for  this 
work  are  two  spectra  recorded  by  Dennis  Keefer  and  his  group  from  the  University  of 
Tennessee  Space  Institute  from  shots  26  and  27  of  the  Thunderbolt  facility  operated  by  the 
Westinghouse  group  headed  by  Hugh  Calvin. 

Thunderbolt  is  an  "experimental  test  bed”  rail  gun.  The  features  that  are  important  for 
this  work  are  shown  in  Fig.  1.  The  inner  part  of  the  Thunderbolt  barrel  consists  of  alternate 
sectors  of  copper  rail  and  boron  nitride  (BN)  insulator.  The  projectile  is  injected  into 
Thunderbolt  at  a  relatively  high  velocity  by  a  light  gas  gun.  When  the  projectile  is  well  into 
the  barrel  a  large  voltage  is  induced  between  the  rails  and  an  arc  is  drawn  through  an 
aluminum  fuse  attached  to  the  back  of  the  projectile.  The  arc  ionizes  the  fuse  and  some  of 
the  gas  thus  creating  a  "plasma  armature".  The  J  X  B  force  on  the  armature  further 
accelerates  the  projectile.  By  the  time  the  projectile  is  several  meters  into  the  barrel,  the 
plasma  armature  contains,  in  addition  to  the  aluminum  from  the  fuse,  a  significant  component 
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of  copper,  boron,  and  nitrogen.  These  constituents  of  the  armature  are  continuously  ablated 
from  the  rails  and  insulators.  Thus  the  set  of  atomic  elements  that  need  to  be  included  in  the 
modeling  are:  aluminum,  boron,  copper,  and  nitrogen.  The  states  included  in  the  atomic  data 
table  are  listed  in  Table  I. 


The  spectra  were  recorded  via  an  optical  fiber  inserted  through  a  radial  hole  in  the 
insulator  (BN).  The  experimental  set  up  is  similar  to  that  described  by  Sedghinasab,  Keefer, 
and  Crowdr*  and  Keefer  and  Crawford^.  The  optical  fiber  used  has  a  band  pass  window  of 
295  -  345  nm.  Although  neither  B  I  or  N I  have  any  strong  lines  in  the  spectral  window 
imposed  by  the  optical  fiber ,  there  are  are  reasons  for  their  inclusion  in  the  atomic  data  base. 
Boron  and  nitrogen  must  be  included  in  order  to  calculate  accurate  opacities.  B I  has  strong 
lines  just  outside  the  window  at  209  and  250  nm  which  are  important  for  energy  transport  and 
wall  loading.  Finally,  there  is  a  weak  B II  line,  2s2p  -  2p^  at  345.1  nm,  just  outside 
the  optical  fiber  window.  The  high  energy  wing  of  this  line  is  in  the  window  and  since  the 
lower  level  of  this  line  is  an  excited  state,  the  opacity  in  the  line  should  be  small  enough  for 
the  wing  to  be  present  in  emission.  The  intensity  and  width  (slope  of  the  wing)  should  be  a 
valuable  diagnostic  of  the  internal  temperature  and  density.  As  an  extension  of  this 
temperature  diagnostic  a  list  of  the  strong  lines  in  the  optical  fiber  window  (OFW)  from 
elements  that  might  be  included  (seeded)  in  the  plasma  armature  is  included  as  Table  II.  In 
this  table  the  first  quantity  is  the  wavelength  of  the  line  in  Angstroms,  the  second  is  the 
relative  strength  of  the  line  in  arbitrary  units,  and  the  third  is  the  emitting  atom  or  ion.  An 
inspection  of  this  table  reveals  several  elements  (Ni  and  V  to  name  two)  that  have  lines  from 
the  neutral  and  first  ionized  ions  that  are  well  removed  form  the  strong  Cu  I  and  A1 1  lines 
and  thus  could  be  valuable  diagnostic  additions. 


COMPUTER  CODE  AND  ATOMIC  DATA 


The  ionic  populations  in  the  plasma  may  be  characterized  by  a  set  of  atomic  rate 
equations  of  the  form 
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where  f-^  is  the  fractional  population  of  atomic  level  i  for  species  m  in  spatial  zone  k,  and 
Wjikm 's  net  rate  describing  the  transition  from  initial  state  j  to  final  state  i.  An  equation 


of  this  type  is  constructed  for  each  of  the  atomic  levels  included  in  the  model.  In  the  present 
instance,  the  assumption  of  collisional-radiative  equilibrium  (CRE)  is  made,  which  involves 
dropping  the  explicit  time  dependence  in  the  above  equation.  The  rate  coefficients  that  are 
used  to  calculate  the  populating  and  depopulating  rates,  are  calculated  using  various 

scattering  techniques  and  the  methods  used  in  calculating  the  corresponding  rate  coefficients 
are  summarized  elsewhere^. 

Radiation  emission  from  the  plasma  and  its  opacity  are  dependent  on  the  local  atomic- 
level  population  densities.  Except  for  optically  thin  plasmas,  however,  the  level  populations 
depend  on  the  radiation  field,  since  optical  pumping  via  photoionization  and  photoexcitation 
can  produce  significant  population  redistribution.  Thus,  the  ionization  and  radiation  transport 
processes  are  strongly  coupled  and  must  be  solved  self -consistently.  In  this  model,  an 
iterative^  procedure  is  used,  where  level  populations  are  calculated  using  the  radiation  field 
from  the  previous  iteration,  then  using  these  populations  to  calculate  a  new  radiation  field  until 
convergence  is  reached. 

The  equation  of  radiative  transport  is  solved  at  a  large  number  of  discrete  frequencies, 
providing  resolution  of  emission  lines,  recombination  edges  and  absorption  edges.  At  each 
frequency  v  and  spatial  zone  k,  contributions  to  the  total  opacity  and  emission  coefficient 
jvk  from  each  bound-b^und,  bound-free  and  free-free  process  are  summed.  For  example,  the 
contributions  from  the  bound-bound  transition  i-i’  of  material  m  to  k^  and  are  given  by 
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respectively.  The  factor  in  brackets  represents  the  reduction  in  opacity  due  to  stimulated 
emission,  nj^  is  the  ion  density  in  zone  k.  E->m  is  the  transition  energy,  g^  is  the  statistical 
weight,  A-’m  is  the  radiative  decay  rate  and  u^^  is  the  line  profile  function.  Voigt  line 
profiles  were  used  in  the  present  calculation.  Similarly,  the  contributions  to  k^  and  from 
the  free-bound  (radiative  recombination)  process  for  atomic  level  i  of  material  m  are 
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for  hv  >  .  The  factor  in  brackets  represents  the  reduction  in  opacity  due  to  stimulated 

recombination.  and  n^  are  the  local  electron  temperature  and  electron  density.  E^  is 
»he  ionization  energy,  is  the  photoionization  cross  section  and  (3q  =  (l/2)h^/(2mnekr'^- 
Finally,  the  contributions  to  ky^  and  jy^  from  free-free  (Bremsstrahlung)  processes  are  given 
by 
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where  Cff=  8ne^/[3mec(6itmek)l^J,  <z^^>  is  the  mean  square  of  the  ionization  and  Gy  is 
the  Gaunt  factor. 

When  the  plasma  is  optically  thick  at  a  particular  frequency  v,  the  emitted  intensity 
will  be  given  jv/ky.  If  the  plasma  is  in  local  thermodynamic  equilibrium  (LTE),  the  intensity 
will  reduce  to  the  Planck  radiation  field  at  that  frequency.  In  LTE,  the  atomic  populations 
satisfy  the  Boltzmann-Saha  relations 
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Then  the  ratios  j^j.  m/k*£  m,  }$£°°fcv  °°,  and  jvfk/I4k 

independently  reduce  to 
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which  is  the  Planck  intensity. 

The  radiative  couplings  Cyj^.  to  each  of  the  other  zones  k’  are  subsequently  calculated 
from  the  opacities  and  the  optical  paths  described  by  the  geometry  of  the  problem.  The  method 
of  obtaining  angle-averaged  escape  probabilities  in  cylindrical  and  spherical  geometries  is 
explained  in  reference^.  The  local  rate  of  energy  density  change  in  zone  k  due  to  radiation 
transport  is  then  given  by 


$v. 


and  the  photoexcitation  rate  for  the  transition  i-i’  is 
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where  k^£  m  is  the  opacity  for  the  bound-bound  process  i-i  defined  above.  Similarly,  the 
photoionization  rate  for  level  i  of  material  m  is 
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where  k^,£  m  is  the  opacity  for  the  bound-free  process  (the  ionization  of  atomic  level  i)  and 
is  the  volume  of  zone  k. 
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For  each  state  included  in  the  modeling,  Table  I  includes  the  energy  level  of  the  state,  its 
statistical  weight,  and  designation  or  name.  The  energy  levels  and  transition  probabilities  of 
ionic  states  were  obtained  from  Cowan’s  code^.  For  atoms  (neutrals),  configuration 
interaction  is  important  enough  that  it  is  more  difficult  to  get  accurate  energy  levels  from 
Cowan’s  code.  On  the  other  hand  some  data  is  in  the  literature,  therefore,  the  energy  levels 
and  transition  probabilities  for  the  states  of  B  I  and  N  I  were  obtained  from  Wiese,  Smith, 
and  Glennon^,  for  the  states  of  A1 1  the  source  was  Wiese,  Smith,  and  Miles**.  The  data  for 
Cu  I  was  obtained  form  references  9-15.  The  collisional  ionization  rates  were  calculated 
using  a  semiclassical  approximation  described  in  reference  16.  The  photoionization  cross 
sections  were  calculated  as  described  in  reference  16  using  the  Karzes  and  Latter*  ^  Gaunt 
factors.  For  dipole  allowed  transitions  the  collisional  excitation  rates  were  calculated  using 
the  semiclassical  approximation  summarized  in  reference  16.  The  distorted  wave  code 
developed  by  M.  Blaha***  was  used  for  the  monopole  and  quadrupole  excitation  rates.  The 
Blaha’s  distorted  wave  with  exchange  code***  was  used  for  those  spin-flip  collisional 
exitations  that  were  included. 


RESULTS  OF  FITTING  OBSERVED  SPECTRA 


The  spectra  recorded  on  Thunderbolt  shots  26  and  27  are  shown  in  Figures  2  and  3.  The 
task  of  the  work  reported  here  is  to  numerically  simulate  these  spectra.  There  are  three 
dominant  features  in  the  experimental  spectra  Starting  a*  3.6  eV  is  the  wing  of  the  B  II  2s2p 
*P  -  2p^  *D  mentioned  above.  Centered  at  3.8  eV  is  a  (partial)  blend  of  two  Cu  I  lines:  4p 
i/2  -  4s  ^>1/2  ^2/2  '  ^l/2-  The  maj°r  feature  is  at  about  4.1  eV  and 

is  a  blend  of  two  A1 1  lines.  These  two  lines  have  a  common  upper  level  and  terminate  on  the 
two  components  of  the  AJ I  ground  state.  These  two  components  are  separated  by  only  0.01 
eV  and  are  treated  as  a  single  state  m  the  model.  Thus  in  the  in  the  synthetic  spectmm  the 
aluminum  feature  is  modeled  as  a  single  line:  3d  -  3p  ^P.  This  latter  line  is  strongly 
quenched  by  electron  collisions  that  excite  the  upper  level  from  the  3d  to  the  4p  level.  The 
3d  and  4p  levels  are  only  0.065  eV  apart  and  the  collision  is  dipole  allowed  so  the  collision 
rate  is  large.  For  Te  =  2  eV,  the  rate  coefficient  is  nearly  10'-*  cm-* /sec!  This  strong 
quenching  results  in  a  very  broad  line.  The  strongest  quenching  of  the  two  Cu  I  lines  results 
from  electron  collisions  that  transfer  the  excitation  back  and  forth  between  the  upper  levels  of 
the  two  lines.  The  upper  levels  of  these  two  lines  are  only  0.03  eV  apart,  but  the  transitions 
are  not  dipole  allowed  and  thus  the  excitation  rate  coefficient  at  Tg  =  2  eV  is  almost  an  order 
of  magnitude  smaller  than  for  the  rate  that  quenched  the  A1 1  line.  Thus  the  width  of  the  Cu  I 
lines  is  much  smaller  than  the  A1 1  line.  These  line  widths  are  demonstrated  in  Fig.  4,  the 
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spectrum  of  an  optically  thin,  homogeneous  plasma  with:  Tg  =  1  eV,  P  =  3  x  10  dynes/cm  , 
and  equal  numbers  of  aluminum,  copper,  boron,  and  nitrogen  atoms  or  ions.  The  dashed 
curve  superimposed  on  Fig  4  is  a  1  eV  black  body  for  comparison.  The  small  dip  in  the 
center  of  the  A1 1  line  profile  is  a  construct  or  artifact  in  the  program.  It  is  the  result  of  an 
interaction  between  frequency  gridding  and  the  overlapping  of  lines,  edges  etc.  Figure  5 
shows  the  spectrum  when  Te  =  2  eV  and  the  plasma  consists  of  .2%  Ai,  20%  Cu,  and  80% 
BN  (%  by  mass).  At  2  eV  Al  I  is  virtually  absent  and  Cu  I  is  reduced  to  the  point  where  the 
wing  of  the  B  II  line  (starting  at  3.6  eV)  now  dominates  the  spectrum  in  the  OFW.  Figures  4 
and  5  show  the  intrinsic  line  shapes  and  widths  at  temperatures  and  densities  characteristic  of 
the  bulk  of  the  plasma  armature.  In  the  remaining  synthetic  spectra  shown  in  this  report  the 
effects  of  opacity,  radiation  transport  and  gradients  in  temperature  and  density  are  combined 
with  these  intrinsic  profiles  to  produce  the  synthetic  spectra. 


A  number  of  calculations  were  made  with  a  multifrequency  radiation  transport  code 
using  aluminum,  copper,  boron,  and  nitrogen  atomic  data.  These  runs  were  made  to  scope 
out  the  composition  of  the  plasma  and  the  temperature  and  density  profiles.  In  these  runs  the 
adjustable  parameters  were:  a)  the  fraction  (by  mass)  of  aluminum  and  copper  (the  remainder 
is  BN),  b)  the  electron  temperature  on  axis  (in  eV),  c)  the  partial  pressure  (dynes/cm^,  of  the 
four  constituents  we  model,  and  d)  a,  the  exponent  in  the  electron  temperature  profile 
equation  suggested  by  Schlichting^  on  the  basis  of  an  analogy  to  the  velocity  profile  in 
rbulent  flow,  and  used  by  Tidman,  Goldstein,  and  Windsor20  for  a  plasma  armature: 

T(ri)  =  Tc(1.0-ri/rw)ct. 

where  rj  is  the  radius  to  the  center  of  the  i  th  cell  and  rw  =  2.8  cm  is  the  radius  of  the  wall. 
Most  of  the  nuts  in  this  set  used  30  radial  cells  and  a  common  radial  grid.  In  order  to  resolve 
the  temperature  gradient  on  the  outside,  the  cell  spacing  decreased  with  radius.  In  Fig.  6  the 
temperature  profiles  for  a  temperature  on  axis  of  2  eV  and  a  =  1/4  (o  o  o),  1/6  (*  *  *),  1/8  (+ 
+  +),  and  1/10  (x  x  x)  are  plotted.  Note  that,  for  the  radial  grid  used,  the  temperature  gradient 
is  well  resolved  and  the  temperature  in  the  30  th  cell  is  small. 


Since  there  are  only  three  spectral  features  plus  the  continuum,  one  might  suspect  that  a 
unique  fit  may  be  "too  much  to  hope  for".  In  addition  the  modeling  is  complicated  by  the 
fact  that  the  experimental  spectra  were  time  integrated,  and  that  the  plasma  armatures  may 
not  be  perfectly  cylindrically  symmetric  and  homogeneous  as  we  assume  in  the  numerical 
modeling.  Although  at  this  stage  we  will  see  how  far  we  can  go  without  them,  time  resolved 
spectroscopy  and  absolute  calibration  of  the  spectrometer  may  be  possible  and  would  ease 
the  modeling  task  considerably. 
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The  opacity  in  the  vicinity  of  the  aluminum  line  is  typically  dominated  by  the 
continuum,  thus  the  width  of  the  aluminum  line  is,  to  first  approximation,  controlled  by  the 
total  mass  density  or  pressure.  With  this  observation  one  can  fix  the  pressure  from  the 
linewidth  at  about  3x10**  dynes/cm^.  Next  one  notes  that  in  the  measured  spectra  the 
copper  absorption  feature  does  not  extend  to  the  base  line.  A  value  of  alpha  of  about  1/10  is 
required  to  accurately  simulate  this  feature  of  the  spectra.  (It  should  be  noted  again  that  we 
assume  a  homogeneous  plasma  with  cylindrical  symmetry  in  the  density  and  temperature,  a 
few  warm  spots  near  the  outside  wall  in  the  actual  plasma  could  easily  lead  us  to  assign  too 
small  a  value  to  a.)  As  indicated,  the  best  fit  of  the  numerical  spectrum  to  that  of 
Thunderbolt  27  was  achieved  with:  Tc  =  2  eV,  P  =  3  x  10®  dynes/cm^,  a  =  1/10, 0.2%  Al, 
20%  Cu,  and  79.8%  BN.  This  spectrum  is  shown  in  Fig.  7  and  8.  Figure  7  is  limited  to  the 
OFW  and  is  a  linear  plot.  Figure  8  is  a  survey  plot  with  the  energy  range  extended  to  1  -  10 
eV  and  the  intensity  is  logarithmic.  Superimposed  on  this  survey  plot  is  the  spectrum 
corresponding  to  a  black-body  at  the  central  temperature,  Tc.  The  survey  plot,  which  clearly 
shows  the  deviation  from  black  body,  is  useful  for  estimates  of  wall  loading  and  radiative 
cooling.  Also  included  as  Fig.  9  is  the  plot  of  the  optical  depth.  The  dashed  curve  on  Fig.  9 
is  the  free-free  contribution  to  the  opacity.  From  Fig.  9  one  should  note  that  while  the  two 
copper  lines  have  optical  depth  of  more  than  1000,  the  wing  of  the  B  II  line  over  600,  even 
the  bound-free  continuum  has  optical  depth  of  100!  Figs.  10  and  1 1  show  the  accessibility  as 
viewed  from  the  outside  (labeled  "Front")  and  from  the  center  (labeled  "Back").  The 
accessibility  is  defined  as  the  fraction  of  the  radius  to  the  point  of  unity  optical  depth.  Thus 
the  outside  ("Front")  accessibility  is  useful  in  determining  the  origin  of  the  radiation  that 
impinges  on  the  wall  and  is  limited  to  the  OFW,  while  the  center  ("Back")  accessibility  is 
useful  in  estimating  the  energy  transported  by  radiation  from  the  hot  interior  outward  and  is 
plotted  over  the  1  -  10  eV  range. 


Inspection  of  the  back  accessibility,  Fig.  11,  reveals  that  even  the  continuum  radiation 
penetrates  at  most  1%  of  the  way  out  of  the  cylinder.  In  view  of  this  fact  one  might  question 
the  sensitivity  of  the  emitted  spectrum  to  the  details  of  the  temperature  profile  at  the  center. 
To  investigate  this  sensitivity  two  runs  were  made  using  the  same  parameters  as  the  previous 
run  but  with  a  spike  (dip)  in  the  temperature  at  the  center.  The  spike  had  a  magnitude  of  2  eV 
and  an  extent  of  0.5  cm.  The  dip  had  a  magnitude  of  1  eV  and  the  same  extent.  The  spectra 
in  the  OFW  of  these  two  runs  were  indistinguishable  from  that  of  the  previous  run  (Fig.  7). 
Thus,  to  nobody’s  great  surprise,  one  concludes  that  the  emitted  spectrum  is  not  sensitive  to 
the  details  of  the  temperature  in  the  innermost  20  -  30  %  of  the  plasma!  Rather  than  proceed 
with  a  set  of  runs  with  temperature  distributions  modified  over  larger  and  larger  portions  of 
the  central  region,  let  us  investigate  the  sensitivity  of  the  emitted  spectra  to  other  parameters. 
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Figure  12  shows  the  spectrum  in  the  OFW  for  the  case  with  Tc  =  3  eV  and  the  other 
parameters  unchanged  from  above.  Note  that  unlike  the  case  described  above,  the  plasma  for 
the  present  case  is  50  %  hotter  everywhere  than  that  for  Fig.  7.  For  Fig.  13  Tc  =  2  eV  but  P  = 
4x10**  dynes/cm2.  For  Fig.  14  Tc  =  3  eV  and  P  =  4x10**  dynes/cm2.  The  sensitivity  of  the 
spectra  to  Tc  and  P  is  evident  from  the  comparison  of  Fig.  7  to  Figs.  12-14.  The  sensitivity 
of  the  spectra  to  a,  the  coefficient  in  the  temperature  distribution  equation  is  evident  from  the 
comparison  of  Fig.  7  to  Figs.  15  and  16.  For  Fig.  7  a  =  1/10,  for  Figs.  15  and  16  a  =  1/8  and 
1/6  respectively.  The  sensitivity  is  evident!  Although  the  figures  are  not  included,  the 
sensitivity  analysis  on  the  composition  has  been  performed.  The  synthetic  spectrum  is 
measurably  different  when  the  composition  of  either  aluminum  or  copper  is  changed  by  25 
per  cent. 


The  radial  gridding  on  all  runs  to  this  point  was  the  same.  In  order  to  resolve  the 
temperature  gradient,  the  cell  size  decreased  with  r.  The  grid  consists  of  30  cells,  with  the 
cell  size  decreasing  as  a  geometric  progression.  The  30-th  cell  goes  from  r  =  2.79991  to  r  = 
2.80000  cm  so  that  the  cell  is  9  10*'  cm  thick.  An  inspection  of  Fig.  10,  the  Front 
Accessibility  shows  that,  at  line  center  of  the  copper  features,  the  accessibility  is  5  loA 
Thus  the  photons  that  strike  the  wall  come  from  a  layer  1 .4  10'~*  cm  thick  and  our  last  cell  is 
optically  thick  at  line  center!  Since  it  was  just  these  line  centers  that  were  used  to  fix  a,  a  run 
was  made  using  a  radial  grid  that  reduced  the  maximum  opacity  through  the  outer  cell  to  1. 
The  other  parameters  for  this  run  were  those  of  Fig.  7.  Comparison  of  this  run  shown  in  Fig. 
17  to  Fig.  7  shows  that  perhaps  we  should  assign  some  uncertainty  to  a.  The  cell  sizes  are  so 
small  that  use  of  the  Schlichting-Tidman  equation  for  temperature  has  to  be  questioned.  One 
may  need  a  wall  loading,  sputtering  etc  model  coupled  to  thermal  conduction  and  (turbulent?) 
hydrodynamic  to  assign  a  temperature  that  close  to  the  wall.  To  follow  that  trend  for  a  bit 
more,  a  couple  of  runs  were  made  in  which  there  is  a  thin  layer  near  the  wall  with  no  copper 
or  aluminum.  The  layer  thickness  was  1  mm  and  0.07  mm  the  resulting  spectra  are  shown  in 
Figs.  18  and  19  respectively.  Figure  18  shows  the  copper  lines  weakly  in  emission,  thus  a 
layer  1  mm  thick  is  ruled  out.  Figure  19,  although  very  different  from  Fig.  7,  can  not  yet  be 
used  to  rule  out  a  0.07  mm  unmixed  layer.  The  question  is:  are  there  parameters  (a,  P,  Tc, 
etc)  that,  when  used  with  that  unmixed  layer  still  yield  a  synthetic  spectrum  that  matches 
Fig.  3?  To  pursue  that  question,  two  more  runs  were  made.  Figure  20  shows  the  effect  of 
changing  a  to  1/6,  and  is  clearly  a  step  in  the  right  direction.  Figure  21  shows  the  additional 
effect  of  reducing  the  pressure  to  1.0  x  10**  dynes/cm2.  Although  the  synthetic  spectrum  in 
Fig.  21  is  now  getting  close  to  that  of  Thunderbolt  shot  27,  that  low  a  pressure  can  be  ruled 
out  on  the  basis  of  other  diagnostics,  thus  one  can  state  that  if  there  is  a  layer  near  the  wall 
where  the  (turbulent?)  mixing  is  not  effective,  it  is  clearly  thinner  than  0.07  mm! 


After  an  extensive  search  of  parameter  space,  none  were  found  such  that  the  numerical 
spectra  matched  that  measured  on  Thunderbolt  shot  26.  The  aluminum  absorption  feature  in 
the  measured  spectra  is  too  narrow  compared  to  the  copper  features.  From  the  above 
negative  results  we  infer  that  the  plasma  armature  on  shot  26  did  not  satisfy  some  of  the 
assumptions,  probably  that  of  homogeneous  mixing. 


CONCLUSIONS 

For  shot  27:  The  spectroscopy  is  consistent  with  the  assumptions  of  cylindrical 
symmetry,  homogeneous  mixing  of  the  atomic  constituents,  pressure  equilibrium,  and  a  very 
strong  temperature/density  gradient  at  the  outside  edge.  For  a  radial  temperature  distribution 
given  by  equation  1,  the  spectrum  is  consistent  with  an  a  of  1/10  - 1/8  and  Tc  =  2  eV. 

Further,  the  spectra  is  consistent  with  P  =  3  x  10^  dynes/cm^,  a  plasma  armature  consisting 
of  20%  copper,  0.2%  aluminum,  and  the  remainder  BN.  If  there  is  an  unmixed  boundary 
layer,  its  thickness  is  smaller  than  0.007  cm.  Although  the  spectral  analysis  can  not  rule  out  a 
central  spike  or  dip  in  the  temperature,  the  requirement  of  homogeneous  mixing  and  strong 
temperature  gradient  on  the  outside  render  it  most  unlikely. 

For  shot  26:  One  capacitor  back  fired  early  on  this  shot,  thus  one  would  not  be 
surprised  if  the  plasma  armature  were  pathological  in  some  way.  From  the  narrowness  of  the 
aluminum  feature  compared  to  those  of  copper,  one  is  lead  to  the  conclusion  that  the 
spectrum  on  shot  26  is  not  consistent  with  some  or  all  of  the  assumptions:  cylindrical 
symmetry,  homogeneous  mixing  of  the  atomic  constituents,  pressure  equilibrium,  and  radial 
temperature  distribution  such  as  equation  1.  Since  the  Cu  I  and  B  II  features  are  similar  to 
those  on  shot  27,  it  is  the  assumption  of  homogeneous  mixing  that  is  most  likely  violated. 

(i.e.  a  "blob"  of  aluminum  near  the  wall). 

Fundamental  Physics:  The  absorption  features  are  sensitive  to  the  temperature  gradient. 
Thus  to  some  extent  one  of  the  quantities  critical  to  the  understanding  of  the  plasma 
armature,  i.e.  the  temperature  gradient,  can  be  inferred  from  spectroscopy.  The  sensitivity  of 
the  spectrum  to  unmixed  boundary  layers  is  an  indication  of  the  strength  of  the  turbulence  in 
the  hydrodynamics. 
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FUTURE  WORK 


Numerical  modeling  of  the  spectra  would  be  enhanced  if:  a)  the  experimental  spectra 
were  time  resolved,  b)  the  instrument  were  absolutely  calibrated,  c)  the  optical  fiber  had  a 
larger  band  pass,  d)  the  plasma  were  seeded  with  elements  that  provide  additional  spectral 
features  (in  the  OFW),  e)  spectra  were  recorded  with  the  optical  fiber  passing  through  the 
rail  rather  than  the  insulator.  The  first  three  are  difficult,  the  last  two  much  less  so.  On  the 
theoretical/numerical  side:  a)  a  better  temperature  distribution,  at  least  very  near  the  wall,  b) 
multi-dimensional  calculations. 
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Table  I 

Atomic  Structure 


ALI 

B  I 

0.0000 

6.0 

3p  2P 

0.0000 

6.0 

2p  2P 

3.1427 

2.0 

4s  2S 

3.3959 

12.0 

2s  2p2  4P 

4.0215 

10.0 

3d  2D 

6.2998 

10.0 

2s  2p2  2D 

4.0867 

6-0 

4p  2P 

7.2436 

18.0 

n=3 

4.7528 

30.0 

3s  3p2 

7.8667 

2.0 

2s  2p2  2S 

4.8270 

10.0 

4d  2D 

8.2794 

32.0 

n=4 

5.1229 

14.0 

4f  2F 

5.8103 

50.0 

n=5 

B  n 

0.0000 

1.0 

2s2  lS 

ALH 

4.6841 

9.0 

2s  2p  3P 

0.0000 

1.0 

3s2  *S 

9.0975 

3.0 

2s  2p  lP 

5.1091 

12.0 

3s  3p 

12.2611 

9.0 

2p2  3P 

11.4262 

4.0 

3s  4s 

12.6877 

5.0 

2p2  *D 

11.4724 

15.0 

3P2 

15.8236 

1.0 

2p2  lS 

12.2336 

20.0 

3s  3d 

183550 

27.0 

2s  n=3  Trip. 

13.1661 

12.0 

3s  4p 

18.8119 

9.0 

2s  n=3  Sing. 

15.2718 

20.0 

3s  4d 

21.8616 

48.0 

2s  n=4  Trip. 

15.3727 

28.0 

3s  4f 

21.9270 

16.0 

2s  n=4  Sing. 

16.4999 

100.0 

3sn=5 

18.2843 

60.0 

3p3d 

B  III 

0.0000 

2.0 

2s  2S 

AL  III 

6.0084 

6.0 

2p  2P 

0.0000 

2.0 

3s  2S 

22.6171 

2.0 

3s  2S 

6.9138 

6.0 

3p  2P 

24.2037 

6.0 

3p  2P 

14.6741 

10.0 

3d  2D 

24.5888 

10.0 

3d  2D 

16.0579 

2.0 

4s  2S 

30.5077 

32.0 

n=4 

18.2215 

6.0 

4p  2P 

33.3282 

50.0 

n=5 

20.9485 

10.0 

4d  2D 

34.8438 

72.0 

n=6 

21.1869 

14.0 

4f  2F 

23.6710 

50.0 

n=5 
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Table  I 
Continued 


Cul 

N  I 

0.0000 

2.0 

*  V 

0.0000 

4.0 

2p3  4S 

1.3890 

6.0 

4*2  2d5/2 

2.2435 

10.0 

2p3  2D 

1.6421 

4.0 

4s2  2°3/2 

3.1185 

6.0 

2p3  2P 

3.7858 

2.0 

4P  2pl/2 

102499 

12.0 

2p23s  4P 

3.8166 

4.0 

4P  2p3/2 

11.6727 

36.0 

2p2  3p  Quart 

4.8376 

6.0 

V%2 

11.8199 

18.0 

2p2  3s  Doub. 

4.9734 

4.0 

V  S/2 

12.8907 

60.0 

2p2  3d  Quart 

5.0763 

2.0 

4P'  V 

132904 

54.0 

2p2  3p  Doub. 

5.1226 

28.0 

4p’  4F 

13.9866 

80.0 

2p2  3d  Doub. 

5.4212 

6.0 

4p'  %2 

5.4710 

20.0 

4p*  4D 

n  n 

5.6780 

24.0 

4p’  Doub. 

0.0000 

9.0 

2p2  3P 

6.1915 

10.0 

4d  2D 

1.9571 

5.0 

2p2  *D 

6.9763 

82.0 

Comp. 

4.8372 

1.0 

2p2  XS 

5.1401 

5.0 

2s2p3  5S 

CuH 

11.4279 

15.0 

2s  2p3  3D 

0.0000 

1.0 

3d10  lS 

133628 

9.0 

2s  2p3  3P 

1.5723 

20.0 

3d9  4s 

18.7515 

9.0 

2p3s  3P 

7.1048 

60.0 

3d94p 

18.9172 

3.0 

2p3s  *P 

12.6703 

100.0 

3d9  4d 

21.1874 

27.0 

2p  3p  Trip. 

15.1924 

140.0 

3d94f 

21.7307 

9.0 

2p  3s  Sing. 

15.6707 

500.0 

3d9  n=5 

23.4969 

45.0 

2p  3d  Trip. 

23.6856 

15.0 

2p  3d  Sing. 

Cu  in 

0.0000 

10.0 

3d9  2D 

N  ffl 

8.3973 

90.0 

3d8  4s 

0.0000 

6.0 

2pl  2P 

15.6580 

270.0 

3d8  4p 

6.9489 

12.0 

2s  2p2  4P 

24.4489 

450.0 

3d8  4d 

12.6293 

10.0 

2s2p2  2D 

29.3734 

630.0 

3d8  4f 

15.8204 

2.0 

2s  2p2  2S 

30.8199 

2250. 

3d8  n=5 

18.0952 

6.0 

2s  2p2  2P 

23.0657 

4.0 

2p3  4S 

33.1757 

18.0 

n=3 

40.8735 

32.0 

n=4 

15 


Table  H 

Strong  Lines  in  the 
Optical  Fiber  Window 


2956.13 

1600 

Til 

3134.11 

2900 

Nil 

3340.34 

1100 

Tin 

2977.13 

1000 

NaQ 

3135.48 

2500 

NaH 

3341.88 

5700 

Til 

2980.63 

1100 

NaH 

3137.86 

1700 

NaH 

3349.04 

4300 

Tin 

2983.570 

1000 

Fel 

3149.28 

2000 

NaH 

3349.41 

2000 

Tin 

2986.00 

1500 

CrI 

3158.16 

6000 

Mol 

335335 

2500 

an 

2986.47 

2100 

CrI 

3162.57 

1000 

Tin 

3354.64 

4100 

Til 

2989.03 

2800 

Bil 

3163.74 

2000 

NaH 

3358.12 

1600 

Mol 

2994.427 

1000 

Fel 

3168.52 

1600 

Tin 

3361.21 

7200 

Tin 

2994.46 

1000 

Nil 

3170.35 

8700 

Mol 

3369.57 

2900 

Nil 

3000.89 

1100 

CrI 

3179.06 

1000 

NaH 

3370.44 

1100 

Til 

3002.49 

4000 

Nil 

3183.41 

3200 

V  I 

3371.45 

4300 

Til 

3003.63 

2200 

Nil 

3183.98 

5300 

V  I 

3372.80 

5700 

Tin 

3012.00 

3700 

Nil 

3185.40 

3800 

V  I 

3377.48 

2900 

Til 

3017.57 

2800 

CrI 

3186.45 

2400 

Til 

3380.28 

1400 

Tin 

3020.639 

1500 

Fel 

3189.79 

1700 

Nan 

3380.57 

3300 
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Fig.  6  Radial  Temperature  Profiles  for  ooo  0=1/4,  * 
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Fig.  8  Theoretical  Spectrum.  For  this  run  run:  P  =  3.0  x  108,  Tc  =  2.0,  a  =1/10,  %  A1  =  0.2,  %  Cu  =  20.,  %  BN  =  79.8 
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Fig.  14  Theoretical  Spectrum.  For  this  run  run: 
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Fig.  15  Theoretical  Spectrum.  For  this  run  run:  P  =  3.0  x  108,  Tc  =  2.0,  a  =1/8,  %  A 1  =  0.2,  %  G 
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Fig.  16  Theoretical  Spectrum.  For  this  run  run: 


/10,  %  A1  =  0.2,  % 
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Fig.  18  Theoretical  Spectrum.  For  this  run  run:  P  = 
with  Cu  and  A1  excluded  from  the  outer  1  mm  shell. 
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Fig.  21  Theoretical 
with  Cu  and  A1  excl; 
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